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Chiral symmetry breaking patterns in the UlIu) x Uji{n) meson model 
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Chiral symmetry breaking patterns are investigated in tlie Ul (n) xUr (n) meson model. It is shown 
that new classes of minima of the effective potential belonging to the center of the Lie-algebra exist 
for arbitrary flavor number n. The true ground state of the system is searched non-perturbatively 
and although multiple local minima of the effective potential may exist, it is argued that strictly a 
UL{n) X UR{n) — > Uv(n) spontaneous symmetry breaking is possible. The reason behind this is the 
existence of a discrete subset of axial-symmetries, which connects various Uv{n) symmetric vacua 
of the theory. The results are in agreement with the Vafa-Witten theorem of QCD, illustrating that 
it remains valid, even without gauge fields, for this certain effective model of the strong interaction. 
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I. INTRODUCTION 



Quantum chromodynamics (QCD) displays an approximate global f/L(3) x ?7fl(3) chiral symmetry. The observed 
'sj" ■ mesonic spectra shows that, this symmetry is broken spontaneously in the ground state of QCD matter. Around 
'the temperature of the scale of QCD {Kqcd), strongly interacting matter undergoes a chiral symmetry restoration, 
■ accompanied by the quark deconfincmcnt transition [23]. Due to difficulties of solving QCD at low energies, effective 
' models built upon the approximate chiral symmetry can give reasonable results regarding the properties of the chiral 
rS |. transition and therefore the phase diagram of the strongly interacting matter. The n-flavored matrix model displaying 
' ' a UL{n) X Uii{n) symmetry [l| give reasonable results for the mesonic spectra and phase diagram IllnlSil. and can be 
extended with quark degrees of freedom (quark-meson model) jHI, Usj and/or the Polyakov loop 16l-l21|. Although 
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(— I . these models cannot be solved exactly either, they provide a simpler theoretical framework than QCD to describe the 
I— '"behavior of the strongly interacting matter at finite temperature and density. The Ua{^) anomaly of QCD can be 
^ I easily included into the model through 't Hooft's determinant term Q, however it is argued that close to the critical 
■ point the instanton density causing the anomaly is suppressed exponentially @. Furthermore, recently new signs of the 
, restoration of the Ua{^) factor around the critical temperature appeared in the literature for n = 2 [2^, therefore close 
^-H ' to the phase transition the anomaly is not expected to play an important role and the pure [/^(n) x Uji{n) symmetric 

' model might have phenomenological relevance. 
CO I It has been known for a long time that there is no infrared stable fixed point of the renormalization group flow of 
] the UL{n) X Ufi{n) model in 4 — e (e > 0) dimensions [2]: e expansion shows that if n > 2, even the existing non-trivial 
T-H .0(2n^) like fixed point becomes unstable in the direction of the quartic coupling. The non-existence of an infrared 
' stable fixed point indirectly shows a transition, which (if exists) is fluctuation induced, and of flrst order. This has 
, been confirmed numerically for the n = 2 case with the use of the Functional Renormalization Group (FRG) approach 
^ ' 0,11, HI], and expected to remain true at n > 2. 

Analysis of the fixed point structure does not tell anything about the symmetry breaking pattern, however, studies 
of the UL(n) x Uu^n) model in the past assumed a priori that the symmetry is broken to Uvin), or in other words, 
^ I the vacuum expectation value of the dynamical matrix field is proportional to the unit matrix. Presumably, the reason 
■ ■ behind this is the Vafa-Witten theorem [3], which states that in vector-like gauge theories with vanishing ^-angles (e.g. 
QCD) and non-zero but equal fermion masses, global vector-symmetries can not be broken spontaneously. Since the 
?7i(n) X Uu{n) symmetric theory is an effective model of QCD, it is plausible to accept that the original statement 
remains true, however, it should not be considered correct without further exploration. This is especially true for 
vanishing quark masses (i.e. no explicit breaking of chiral symmetry), for which the theorem only states that if a 
vector-symmetry breaking state existed, its energy could not be lower than the symmetry conserving one's. The 
situation is getting more puzzling, if we accept the results of (23 |. where (in the large-n limit) a new local minimum of 
the energy functional realizing a symmetry breaking pattern UL{n) x Unin) — > Uvin — 1) was found, from which the 
system can undergo a first order transition with appropriately chosen model parameters (25j . Furthemore, in (26l. \2l\ 
it was conjectured that Yukawa couplings of the lepton masses of the Standard Model might be described by the 
spontaneous breaking of a hypotethic SU (3) x SU (3) symmetry, realizing yet another ground state which breaks the 
symmetry to SU{2) x SU{2) x U{1). This solution has already been reported in [EJ for a generic number of flavors. The 
existence of these different vacua in the UL{n) x Uii{n) model might be in apparent contradiction with the Vafa-Witten 
theorem, therefore it is an open question whether these minima are not physically applicable for some reason, or the 
applied approximations are inconsistent. 
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The goal of this paper is to answer the aforementioned questions and to map the vacuum structure in an extended 
space of the condensate space, searching for new local minima of the effective potential for arbitrary flavor numbers. 
The paper is organized in a pedagogical way and goes from the most simple cases to a general analysis. In Sec. II, we 
introduce the model and the symmetry breaking we would like to study in the first place. We list the group invariants 
and discuss the possibility of breaking of the diagonal Uv{n) symmetry. In Sec. Ill, we solve the field equations at 
classical level, and show the energy properties with the corresponding spectra of the found solutions. In Sec. IV we 
include quantum fluctuations using Functional Renormalization Group (FRG) formalism. We start with the n = 2 case 
in the so-called Local Potential Approximation (LPA), then generalize the arguments for the exact FRG formulation 
for n > 2. It is followed by an analysis of having a most general condensate, and a discussion of obtaining all the J7y(n) 
symmetric ground states. Conclusions can be found in Sec. V. 



II. THE MODEL AND SPONTANEOUS SYMMETRY BREAKING 



The dynamical variable of the J7l("-) x Un {n) scalar theory is a matrix field M, as an element of the U (n) Lie-algebra. 
It can be parametrized as: 

M== (s'^ + i7r'')r", (a = 0...,n2 - 1), (1) 

where s" and tt" are the scalar and pseudoscalar fields, respectively, with T" being the group generators satisfying 
Tr {T'^T'') = (5"V2- The Lagrangian of the model is 

C = d^Md'^M^ - Tr {MM'^) - ^[Tr {MM'^)f - ^ Tr (MM^MMt). (2) 

In order to have a classical potential bounded from below, we require g\ + .92 > and 52 > [3]. It is clear that, ([2]) is 
invariant under the symmetry transformation 

M UrMUI (3) 

Ul and Ur are independent U{n) matrices realizing the chiral symmetry with parameters and These transfor- 
mations can be reformulated in terms of vector- and axialvector-symmetries as: 

M — > VMV\ M — > A'^MA\ (4) 

where V and A are also unitary matrices with parameters Oy = {9]^ ± d1)/2, respectively. Even though in QCD the 
axial symmetry is broken anomalously at quantum level, in this paper we do not investigate quantitatively the role 
of this anomaly, nor any explicit symmetry breaking. Following [24| . first we are interested in the possibility of the 
symmetry breaking pattern 

<M>^ t7oT° -I- vsT^, (5) 

which means that we set all the expectation values in advance to zero, and investigate only on the vq — wg plane. We 
will come to the question of general condensates later. In ([5]) is the generator of the U{1) subgroup of U{n) and 
refers to the longest diagonal generator, regardless of the actual number of flavors (see also Appendix A). This notation 
is just an adaption from the phenomenologically most important n = 3 case. Note that, pattern ([5]) does not apply 
for n = 2, since in this case there are only four generators. Besides T°, only of them is diagonal: T^, the third Pauli 
matrix. Formulae given below can be applied to n = 2 by considering instead of T®, i.e. substitute index '3' instead 
of '8' everywhere. 

In terms of the coefficient fields and 7r°, ([2]) turns into 

£ ^ ia^s-a'^s- + ^d^TT^dV - V{s^,7r-), (6) 



where 



with the use of notations 



V = m'h{s-,n-) + + -/2(.s^7r'^), (7) 

n 



/i(s°,7r'') = Tr(MAft) = i(s'^s"-f tt^tt'^), (8a) 
his'^n'') = Tr(MMtMMt) = id'''''=d''"^(s^s'= + 7rV)(s^s-'^-K7r%-^) + ir''"r'=^s''s%%^. (8b) 
Here d'^'"^ and /'^'"^ are the totally symmetric and antisymmetric tensors of the U{n) group, respectively, see details in 
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Appendix A. We will need the invariants and their derivatives in a background of vo,vs. The invariants are: 



1 



/2n 



1888 3 

i VsVq. 



(9) 



The needed structure constants are listed in Appendix A, the derivatives of the invariants can be found in Appendix 
B. The mass matrices are defined as follows: 



'3 ds'dsi 



M. 



2(7.) 



(10) 



At classical level these can be written as 



2{s) 



M, 



2M 



m 



251 
251 



dsidsj 



diTidni 



92 92/2 



n dsidsj 
92 d^h 



2gidh 

•n? dsj 



dh 
dsi 



(11a) 
(lib) 



Checking the detailed expressions of (|lla|) and (lllbp found in Appendix B, it can be seen that we have to diagonalize in 
the — 8 indices. This leads to the multiplet structure 1 ® 1 ®n{n — 2) © 2(n — 1) in each sectors. Note that, when n = 3 
these groups correspond to 77, 77', three pions and four kaons in the pseudoscalar sector, and a, fo singlets with three Oq's 
and four k's in the scalar sector. In general, the spectrum may show the possibility of a [/^(n) x Un{n) — > Uv{n — 1) 
symmetry breaking. We need to investigate if this indeed realizes in the minimum of the effective potential. 



III. NEW CLASSICAL SOLUTIONS OF THE FIELD EQUATIONS 

In this section we present the field equations for the condensates v^^Vg, and search for all possible solutions of them. 
We note that the formulae below can be applied to n = 2 as well, we just have to consider index '3' instead of '8' 
everywhere. Using (|B1[) . the field equations are: 



= 



= 



dV 
dV 



\ — 1 

,2„ 



nv 



! , .91/ 2 , 2n , 92 („2 o 2 n^-3n + 3 3\ 
v% + -^(vq + Ug)w8 + 3woW8 - 3 v^vo H ■ (12b) 



Being cubic, these coupled equations can be solved analytically and we obtain nine solutions. Since the potential reflects 
(wojUs) < — > {~vaT—v^) symmetry, we have only five inequivalent solutions. Not considering the parity pairs, we get 
(see also Fig. 1): 



- 0, = 0, (13a) 







[//] _ -nrm^ [//] 



„["fl _ o^./--^ 4^^^1^2j'"^'(""^), (13c) 




0, (13b) 



91 +92 



[IV] _ I -m^n jiv] _ I -m^n{n-l) 



9i{n - 1) + 92n' 



(13d) 



jy] _ / ^[v] _ _ -m'^{n - 1) 

9i/n + g2' * y 9i/n + g2 



We can immediately see that, for n — 2, solutions (|13p are symmetric under the interchange of the two condensates. 
This is due to the fact that at n = 2, there is a vq i — > vs (or more appropriately vq i — > W3) symmetry of the Lagrangian 
and of the field equations as well: e.g. solution [V] is the refiection of solution [IV] with respect to the vs = line (this 
is not true for n > 2). 

Solution [V] shows <M>~ diag (0, 0, 1), which has already been reported in [5j and particularly for n = 3 in 
[26l[27|. where it has been conjectured to be a ground state which may describe lepton mass hierarchy in the Standard 
Model. We would like to mention however that, in the stability region (i.e. 91 + 92 > 0, 92 > 0) this solution always 
describes a saddle point (this can be verified by checking the mass matrices |B3| ). and therefore not a local minimum of 
the potential. In general, solutions [II] and [III] are local minima, [I] is a local maximum, the others are saddle points. 
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3i=1, 92=1 
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FIG. 1: Solutions of field equations p2p for arbitrary n > 2. As n increases, the solutions (except [I]) get away from the origin. 
Solution [I] is a local maximum, [II] and [III] are local minima, [IV] and [V] are saddle points. 



We note that, if the couplings satisfied 52 < and gi + ng2 > 0, the only local minimuni would be solution [V], but 
since in this case the potential is not bounded from below, we do not investigate it further. 

In the large-n limit, solutions [II]-[III]-[IV] (after rescalings Vq — >■ nvo, vs — >■ \/nvs) exactly reproduce the three local 
extrema found in [iSj- In this case solutions of vq become degenerate, and a picture of a double- well potential formes 
in the > direction. The appearing new local minimum (solution [III]) was a candidate in [25'| to describe the 
ground state and a first order transition of the system, see also Fig. 2. Our general analysis shows that this structure 
is although distorted at finite n, it still exists: a new local minimum (solution [III]) emerges for every each n value. 

The energies of the solutions are as follows: 



vii'K.i'^ = 0' n,i"^ .1"! = -77— y\ 



4(51 +.g2)' %i+52)' 
m^ri^(n— 1) rn^n? 

[IV] = -— ; r , ^U^l Av] = -77 ; T, (14) 

and the following relations hold between them (using obvious shorthand notations): 

V[ii] = V[///] < V[iv] < V[v] < V[iy (15) 

Equality between V[iv] and Vjyj is only possible for n = 2, as it should. We see that having two minima with exactly 
the same energy is a general property for arbitrary n (it is interesting to mention that, the energy over the number of 
degrees of freedom is even n- independent in these minima). This could make us suspicious, whether there is a symmetry 
forcing these minima to have the same amount of energy. Let us check the spectrum in these minima as well. First we 
calculate the masses which belong to [II] : 

<[//]- 0, 0, M|^[,,] = -m^-^, (16a) 

= Km=-^'^^^^ ^ ^"''^r?^' (16b) 

and then the ones corresponding to [III]: 

Mliiii] = 0, Af2 = 0, Ml^jjj^ = 0, M|,^[,,,] = 0, (17a) 

It is interesting to note that only excitation multiplicities depend on n, masses do not. We see that ([T7| describes the 
usual UL{n) x UR{n) — > Uv{n) breaking, but the spectrum belonging to p6)) is the same! We just have to interchange 
K i — > K, but apart from parity properties that is just a matter of naming the excitations. Since the energies of 
these states, and also the spectra are the same, both describe the very same symmetry breaking, there is no physical 
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1, 92=1 




FIG. 2: Potential energy at n = oo as a function of vs after rescaling vs — > y/nvs- The value of vo is set to where local minima 
were found. The energy is shifted by a constant to take its zero at vs = 0. Surprisingly, the two minima describes the same 
symmetry breaking pattern. 

distinction between them. Therefore, we can conclude that, at least at classical level we found two physically equivalent 
states, and it cannot be distinguished, which one is the ground state. At this point, we would like to emphasize that, in 
the large- rt limit it can be very misleading that, extrema [II] -[III] -[IV] lie on the same line of vq/ti = \J —"rn^ I {g\ + 52) 
and while minimum [II] has — 0, [HI] shows 7^ (Fig. 2). Since the value of vq is the same for [II] and [III], one 
would expect that obviously these correspond to different symmetry breaking patterns, but eventually it turns out not 
to be true. 

IV. QUANTUM LEVEL 

We would like to make connections with previous works, therefore we choose to include fluctuations in the framework 
of the Functional Renormalization Group (FRG) formalism 6]. In FRG, one derives a flow equation for a one parameter 
(fc) quantum effective action F^, which interpolates between the classical action (describing micro-physics) and the 
quantum effective action including fluctuations. For some general bosonic field(s) the one parameter Schwinger 
functional in euclidean space is defined as 



-Wfc[.7] _ 



(18) 



where J is a source and Rk is a regulator function suppressing fluctuations with momenta lower than fc. The fc dependent 
quantum effective action is defined as: 



rfe[$] = Wfc[j] - / J • $ - i / $ • • 



(19) 



which obeys the so-called flow equation: 



^ ^ i Tr 
dk 2 



(20) 



(2) 

where the trace has to be taken in both functional and matrix sense, and F^ is the second functional derivative of F^. 
The flow starts from a UV cutoff scale A, where by definition (5 being the classical action) 

Ffe=A = S, (21) 



and goes to zero including all the fluctuations, eventually reaching the quantum effective action F^^g = F. 

A common approximation scheme is the Local Potential Approximation (LPA) , where it is assumed that F^ has the 



following form: 



= f {d,^x)dMx) + Uk{x- $)) , (22) 



where Uk is the fc-dependent effective potential. This is expected to be a good approximation if the anomalous 
dimension remains small. With the choice of the commonly used regulator function introduced by Litim : Rk (q) = 
(fc^ — q^)Q{k'^ — q^), from (PO)) . in 3+1 dimensions we get 

dUk ^ 1 



dk 67r2 k"^ + 

i 

where Mi are the eigenvalues of the mass matrix d^Uk/d^id^ 



XI i.2 I A^2 ' (23) 



A. Analysis of the n = 2 case 



We start to analyze the quantum effects in LPA at n = 2, which is the simplest case. Since the U{2) group has only 
two independent invariant tensors (i.e. the previously introduced /i and the ^-dependent effective potential is a 
two-variable function: Uk = Uk{Ii, 12)- Using Appendix A, the mass matrices can be easily built up. In general we 
have: 



2is) 



dlUk 



M, 



+ d^Uk^^ + dnUk — — + d22Uk — — 



di2Uk 



dhdh^ dh dh 

dn^ dn^ dn^ dn^ 



(24a) 
(24b) 



With the help of (jBl[) and (|B2I) . we can evaluate these matrices at (vqjVs), and after diagonalization in the — 3 
sectors, we recover the results of 28]: the right hand side of indeed depends only on Ii and h, as it should. The 
most important observation is that h and h are invariant under the transformation vq -s — > W3, therefore so does the 
right hand side of (1^^ . If we start to flow the effective potential, this leads to the observation that since Uk=A itself 
has vo < — > V3 symmetry, it is preserved throughout the flow to fc = 0. In other words, on the two-dimensional plane 
(vqjV^), the line V3 — can be exactly mapped onto the line vq — 0, which leads us to the fact that, the two local 
minima of (jl3l) evolve along these lines exactly the same way, being undistinguishable physically even at the very end 
of the flow (i.e. the quantum effective action). This concludes that even though we observed a potential symmetry 
breaking scenario Uhin) x Un{n) — > Uyiji — 1), at least at n = 2 this turns out to be UL{n) x UR{n) — > Uv{n) even 
at quantum level. Note that, in 0, H, [13] it is shown numerically that, the system undergoes a first order transition 
starting from the minimum of the line 1)3 = 0, therefore it does so from the minimum of the line uq = 0. 

Since we used LPA, and restricted ourselves to the n = 2 case, we have not yet reached our final goal: we need to 
analyze the situation carefully for n > 2 without having any approximation. The previous analysis however can easily 
be generalized, leading us to an exact result. 



B. Analysis of the n > 2 case 



We start from the exact fiow equation (pO)) . The right hand side is U{n) x U{n) invariant by construction, but this 
time we have more than two invariants, which will be introduced by fluctuations. Since we are interested in the effective 
potential at homogeneous field configurations, we do not have to care about invariant terms containing field derivatives. 
The terms which are to appear and relevant to our interest in the effective action are the following traces: 



Tr(MMl'), Tr(AfMtMMt), Tr (MM^MM+MMt), 



(25) 



Depending on the value n, we have various numbers of independent invariants evolving throughout the fiow. The 
important observation is that all of these are made out of traces of products of M^M. Let us take a closer look at this 
entity in a background uo,W8- Since generators and are real, M|^Q_„g = M^l^^^^g, therefore 



'0 1 



1/ 



1 



2n{n - 1) 



'0 1 



i diag [ (.0 + ^) , (.0 + , (.0 - V^lvsY 



V8 



-{n-l)/ 



(26) 



(27) 
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The important observation here is that, there are two choices to obtain MM^|„Q.t,g ~ 1, which obviously describes a 
symmetry breaking [/^(n) x Uji{n) — > Uv{n): 



vo + = ±(i;o - V^^vs). (28) 

yn — 1 

The plus sign leads to — 0, which is the obvious choice, but the non-trivial scenario is obtained choosing the minus 
sign: 



vs = 2^ -vo. (29) 

n — 2 

(For n = 2 we get vq — 0, as we should.) Checking (|13cp (solution [III]), it lies exactly on this line. The existence of 
that minimum is therefore a consequence of this symmetry, which will be formulated more precisely in the next section. 
We see that for every point of the — line can be associated with a point of the line vs = in the sense that 

MM^ does not change. Because of the right hand side of ([20|) has to depend only on MM^ (and not M itself), the 
aforementioned two lines evolve exactly the same during the flow. Furthermore, because of continuity, not just the lines 
themselves, but also an open neighborhood of them have this property. This leads us to the observation that, along 
these lines not just the value of the potential, but also the appropriately chosen directional derivatives are the same. 
That is what we saw before explicitly; at classical level the calculated spectrum was the same, just some directions of 
the (s°,7r°) space had to be interchanged. This result shows that the applied heavy scalar approximation of (23. [25j 
might look reasonable, but obviously hinders the way in which the exact result was obtained: without having the scalar 
excitations included, the UL{n) x UR{n) invariance is broken in (pO)) . 



C. Minima of the effective potential for a general condensate 

We can formulate the appearance of the new minimum and the symmetry transformation behind it in a more general 
way, and argue that, it is a discrete subset of axial symmetries which generates new minima. Let us start from a 
condensate proportional to the unit matrix: <M>= vqT^ . Since T° commutes with everything, it is clear that Uv{n) 
transformations leave <M> invariant, while UA{n) does not: 

<M> — > V <M> V'l =<M>, <M> — > At <M> A'' ^<M> {A'if. (30) 

The effective potential is symmetric under Uv{n) x C/4(n) rotations, but since vector symmetries leave <M> invariant, 
we need to search for appropriate axial-symmetry tranformations to get copies of the Uv (n) symmetric vacuum. With 
the notation U = {A'')'^, the most general transformation we can apply is <M> — ><AI> U, where U is an arbitrary 
unitary matrix. This immediately shows that not just the invariant tensors, but <M><Mt> itself has to remain 
unchanged, in agreement with the previous subsection. We would like to search for copies of the <M> background, 
where the condensate lies in the center of the Lie- algebra. This means we search for U matrices for which exist {va} 
real numbers fulfilling the following equation: 

<M>U^ ^aT". (31) 

a— diag 

This can be rewritten as 

a— diag 

where aa = sphiVajvo. This shows that, the intersection of the center of the Lie-algebra and the U[n) group itself in the 
fundamental representation connects various non-trivial vacua of the effective potential: each of these transformations 
refer to a copy of the original \Jy(n) symmetric minimum. Unitarity implies: 

1 = f/[/t = Y aaabT^TK (33) 

a.b— diag 

Since T'^T'' = 2(d"''= + ip^'')T'', using that f^" is totally antisymmetric and r° = ^1, jSa]) is equivalent to the 
conditions 

f X! ^ = X! OiaabcP^", (34a) 

a— diag a, 6— diag 

which come from linear independency of the generators. looks simple, but from a practical point of view, it is 

better to express unitarity in terms of matrix elements of UU^ . Using the definition of diagonal generators given in 
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Appendix A, from (l33l) we get: 

n-l 



diag 



'2n 



E 



ao 



n-j-l 



E 



(35) 



We just showed the zeroth and the j*'' element of the diagonal [j E {0...,n — 1}). Requiring each element to be equal 
to unity, we get n equations for n variables {aj}. 

We present a possible method to solve these coupled equations. The {j — 1)*'* and the j*'' equation look as 



ao 



/2n 



V^jU - 1) 



a 



j+i 



y — 

n-j-l 



E 



If j 7^ 0, we can plug (j36bp to ([36a|), and arrive at 



= 1 



(36a) 
(36b) 



a 4 



2j 



2j 



± 1 = ±1, 



(jG {!..., (n-l)}), 



(37) 



which is a recursion for aj and where the ± sign on the left hand side is not arbitrary, it depends on the sign of (j36b[) . 
i.e. the one in the equation. Because of that, it is convenient to start the recursion backwards. We have not yet 
used the last matrix element of (|34p . Requiring it to be equal to unity, we get the starting point of the recursion: 



2n 



y/2n 



\Jn - 1 = ±1. 



(38) 



In practice we start with ((38|). choose a sign and express a„_i with ao. Then we plug the expression of a„_i into ([37]) 
(sign on the left hand side has already been determined with the previous choice!) for j = n — 1, and choose a sign 
again for the right hand side. The procedure continues until j = 1, when a\ is determined explicitly, and with that, 
we get the expression of ao and therefore all the {a^} immediately. 

The solution described in the previous subsection refers to a set of choice, where at each but the first step -plus sign 
is chosen: it leads to 



vo = vo[n-2)ln, Vn-i ^ 2voVri~^/n, Vi = Q (z = 1..., n - 2), (39) 

which is consistent with (|29l) or (jl3bp with (I13cp . (Note that previously index notation "8" was used to [n—l] regardless 
of the value of n.) Choosing plus signs only leads to the trivial 

vq^vq, -D, = (i = 1) (40) 

solution. 

With the procedure outlined above, since there are 2" ways of choosing signs, we can find in general all the 2" 
solution of ([55t very easily, which connects the usual vq ^ Q,Vi — Q line with appropriately chosen directions, and 
which give every local minima of the potential ([7]) on a diagonal background. By construction, these turn out to be 
Uv{n) symmetric ground states, and more importantly, evolve identically with the FRG fiow. This illustrates that 
vacua appearing at classical level can not break the Uv{n) symmetry, and because of the discussed invariance of MM\ 
this property remains true after including quantum fiuctuations too. We note however, although it is not expected, but 
cannot be excluded that, higher order operators may develop even more new local minima not appearing at classical 
level. 

Finally, let us make a comparison between the Vafa-Witten theorem Q and our results. The theorem claims that, 
in vector- type gauge theories with zero 0- vacuum angles (e.g. QCD), if the fermion masses are equal and non-zero, 
none of the vector symmetries can be broken spontaneously. Our result mainly agrees with this, remarkably without 
the appearance of any gauge fields, in an effective model. Adding the following explicit symmetry breaking term to the 
Lagrangian [H = haT'^/A): 

Ch - Tr {H{M + A/t)) = has\ (41) 

paper [1] proves that if H ^ 1 , which refer to equal quark masses, pattern of the symmetry breaking is Ul (n) x Ur (n) — > 
Uvin), moreover, it claims that, if _ff = and because of some reason new, maybe even vector-type symmetry breaking 
minima existed, their energy cannot be lower than the energy of the symmetry conserving state. Our analysis extends 
this argument and excludes the possibility of having a Uv{n) breaking state. In the chiral limit we do have multiple 
ground states as ^4;] suggests it as a possibility, but always the pattern UlIti) x Ua^n) — > Uv{n) occurs. However, 
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in general, explicit symmetry breaking may split these states, therefore H can create physically inequivalent 
stable and/or metastable states. For example ~ 1 with ho > would clearly select solution [II] (jl3bp as the true 
ground state. Phenomenology shows IT] that for QCD ho > 0, hs < 0, with ss 0, which makes (jl3b|) even more 
favorable energetically. Treating the model as an effective theory of strong interaction, this points out the correctness 
of investigating the <M>— wqT" minimum in the first place, but it is important to stress here is that the vacuum 
structure of the model shows a much richer structure than it has been investigated before. 



V. CONCLUSIONS 



We investigated the possibility of having a spontaneous symmetry breaking pattern of the UL{n) x C//j(n) meson 
model other than Uvin). The Vafa-Witten theorem [4] states that for vector- type gauge theories with vanishing 9- 
vacuum angles and nonzero but equal fermion-masses, none of the vector-like global symmetries can be broken. Since 
this applies to QCD and the UL{n) x Uii{n) model is an effective theory of the strong interaction, it is common to 
consider that the expectation value of the dynamical matrix field is <M>~ 1. However, we have no reason to believe 
that, the model obeys a theorem valid for QCD, without further investigations. 

Following the route of [13, HH], we considered a more general, <M>= vqT^ + v^T^ type of symmetry breaking, 
and started our study at classical level. We found that new syrnmetry breaking minima of the effective potential 
appears for arbtirary n, also reproducing the large-n result of [31 and a solution given in [1, [H, US]- Even though 
the examined pattern is capable of describing the breaking of vector-symmetries, it turned out that, in the observed 
minima of the effective potential, all of them were conserved. Furthermore, all the states have the very same energy 
and spectrum, therefore are physically equivalent to the "usual" minima of the effective potential described by the 
condensate <M>~ 1. 

We checked if the structure remained the same after including quantum fluctuations, with the use of the Functional 
Renormalization Group approach. We started with an approximate solution, the Local Potential Approximation for 
n = 2, and found that fluctuations did not change the structure had been described before. Then, generalizing the 
applied arguments, we showed that the results hold even without any approximations, for arbitrary flavor number n. 
The reason behind this turned out to be the existence of a discrete subset of the axial-symmetry, which connects several 
Uv{n) symmetric vacua. This subset can be formulated as the intersection of the center of the U{n) Lie algebra and the 
group itself. Using this symmetry, we managed to show how to obtain all the possible vacua of the theory considering 
the most general diagonal condensate. These minima of the effective potential are Uyin) symmetric by construction, 
which illustrates the validity of the Vafa-Witten theorem for the UL{n) x Unin) meson model. 
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Appendix A. PROPERTIES OF THE U(N) GROUP 

The U{n) group has n independent diagonal and n{n — 1) independent non-diagonal generators leading to in 
total. They are traceless and normalized as Tr (TT^) — S°'''/2. We introduce the generalized Gell-Mann matrices as 
A* = 2r*, where T* are the group generators. The diagonal ones read as: 



A(i 



1 



/I 



V 



.../ 



n{n — 1) 



/I 



V 



-{n-l)) 



(Al) 



We note that, except for the last subsection of section IV, throughout the paper conventionally "8" index goes to 
^("-1)^ independently of the actual value of the flavor number n: 



2T8 = A« = A("-i). 



(A2) 



The non-diagonal generators form two groups, as generalizations of {x) and (y) type of Pauli matrices in the following 
sense: 



bk- 



(A3) 
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Here compact index notations {x,jk) and {y,jk) {j < k) were introduced. The structure constants of the algebra are 
defined through the relation 



A"A^ = {dabc + ^/afcc)A^ 

where dabc is totally symmetric and fate is totally antisymmetric in their indices. We have: 

{A^A^} = 2dab,A^ =^ Tr[{A^A''}A=] =4d,6„ 
where [., .] and {., .} refer to commutation and anticonimutation, respectively. Alternatively, we can write 

fabc = i3Tr(A"A''A^), dabc = i5RTr(A''A''A^). 



(A4) 

(A5a) 
(A5b) 

(A6) 



The structure constants needed for obtaining the mass matrices are: 



Mab 



n 



j888 



(2-n) 



(2 - 2nin~i) ^'^^ if ^ ii^^jn), {y,jn)} 



i(n-l) 



n 



n{n — 1) ' 



else 



(A7) 



Appendix B. CLASSICAL MASS MATRICES 



In this section we show the exact form of the mass matrices at classical level. First we list the first and second 
derivatives of the invariants on a background vq, needed for the expressions of the masses. The first derivatives: 



dh 



0, 



dh 
: 0. 



-vlS^'' + -^d^^\lv^ + -vlv^d^"" + + -vl5''\ (Bla) 

n n 1 n 

(Bib) 



The non-zero second derivatives: 



d'^h 



ds^dsi 



d^h 



dir'-diT^ 



(B2a) 



ds^dsi 



d^h 



-{3vl + vl)6'^ + sJ-d^'^VoVs + id888^8.,„2 _^ l^S'SpS ^ ^.0^.8)^2 ^ (^8«)2^^,#0,8^2^ 



+d' 



(B2b) 
(B2c) 



Relations (|Bip are also needed for the derivation of the field equations, (|B2p can only be used for determining the 
mass matrices. Using the expressions of the structure constants (jA7p . invariants (0) and their derivatives (jBl[) . (jB2p 



dir'-diT^ 



I'D ■'"8 



n 1 \ n ^ — ' 
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just listed, the nonzero elements of the mass matrices in a background vq, are as follows: 



m 



'm' + ^ {gi{ivl + vl) + ig2{vl + 

5iK'+3«|) + 3,g2(«o' + - 
9M+vl)+92{H + - 
.91 ("o + 4) + 3.92(^0 +4 + 7=t""^8)) , 
m' + ^ (51 + 92) {vl + vl) , 



-3n+3 „,2 
: — 1 t^s 



-3k+3 „,2 
: — 1 



3^ 



fool's. 



TO 



'5i(wo+«i)+ 52(^0 
> 51 (wo +52(^0 
i5i(wo +«i)+52(wo 



-3n+3„,2 



n-1 "8 
'-K+l„,2 



n-2 



n-1 



i = 0, j 
i = 0, j = 8 > 
1 = 8, j = 8 
i=je {{x,jn),{y,jn)} 
i = j, else . 

1 = 0, j=0 

i = Q, j = 8 or 1 = 8 

1 = 8, j = 8 

i= j e {{x,jn), {y,jn)} 

i = j, else . 



i = 8, j=0 



(B3a) 



J=0 



(B3b) 
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